NASA Technical Memorandum 88368 


/M -3^ 


Airborne Astronomy Program 
Medium Altitude Missions Branch 
Preprint Series 063 


(KASA-lll-8836e) THE ICMZAllCK SlfiOCTOEE Of N87-19313 

lEE OBICH KEBCLA: JNFBABED IlKi C ESEfiVATICHS 
AfcC (3CDELS (HAS A) 46 p CSCL 03A 

[Jnclas 

G3/89 43769 


The Ionization Structure of 
the Orion Nebuia: Infrared 
Line Obsenfations and 
Models 

J.P. Simpson, R.H. Rubin, 

E.F. Erickson, and , 

M.R. Haas \ 


October 1986 



IVIASA 

National Aeronautics and 
Space Administration 



NASA Technical Memorandum 88368 


Airborne Astronomy Program 
Medium Altitude Missions Branch 
Preprint Series 063 


The Ionization Structure of the 
Orion Nebula: Infrared Line 
Observations and Models 

J. P. Simpson, Lick Observatory, University of California, Santa Cruz, Califronia 
R. H. Rubin, University of California, Los Angeles, California 
E. F. Erickson, Ames Research Center, Moffett Field, California 
M. R. Haas, Mycol, Inc., Sunnyvale, California 


October 1986 


ru/\sA\ 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Field, California 94035 




ABSTRACT 


Observations of the [0 III] 52 and 88 um lines and the [N III] 57wm line 

have been made at 6 positions and the [Ne III] 36 ym line at 4 positions in the 

Orion Nebula to probe its ionization structure. The measurements, made with a 

-40” diameter beam, were spaced every 45" in a line south from and including the 

Trapezium. The wavelength of the [Ne III] line was measured to be 

36 . 01 3±0.004ym. Electron densities and abundance ratios of N**/ 0 *^ have been 

calculated and compared to other radio and optical observations. Detailed one 

component and two component (bar plus halo) spherical models were calculated for 

exciting stars with effective temperatures of 37~40,000K and log g -4.0 and 4.5. 

Both the new infrared observations and the visible line measurements of oxygen 

and nitrogen require T^^^ < 37.000K. However, the doubly ionized neon requires 

a model with T^^^ 2 39,000K, which is more consistent with that inferred from 

the radio flux or spectral type. These differences in T^^^ are not due to 

effects of dust on the stellar radiation field, but are probably due to 

Inaccuracies in the assumed stellar spectrum. The observed N**/0** ratio is 

almost twice the N*/0* ratio. Our beat fit models give N/H-8.4x1 0~^, 0/H-4.0x10”^, 
-4 

and Ne/H-1.3x10 . Thus neon and nitrogen are approximately solar, but oxygen 

is half solar in abundance. From the infrared 0** lines we conclude that the 
ionization bar results from an increase in column depth rather than from a local 
density enhancement. 

Subject Headings ; Nebulae: H II Regions - Nebulae: Orion Nebula - 
Nebulae : abundances 
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I. INTRODUCTION 


H II regions are studied in order to learn about interstellar elemental 
abundances and the far ultraviolet spectra of hot stars. The Orion Nebula is a 
case to study because it is nearby, bright, reasonably compact, the extinction 
is not large in the visible and is negligible in the infrared. There are radio 
maps available at high spatial resolution for determining the electron density. 

It has long been studied in the radio and visible, but only recently in the UV 
and infrared. 

Infrared forbidden lines are important in an observational study of the 
structure of an H II region because some of them arise from ionization states 
that do not have optical forbidden lines. Such lines include transitions of N*^ 
at 57uffl« C* at 158ym, Ne* at 12.8pm, at 10.5ym, and Ar* at 6.98ym. 

Additional advantages are that the lines are less affected by extinction than 
are the optical lines and not at all by scattering. Since infrared forbidden 
lines (unlike optical) are Insensitive to electron temperature, the abundances 
determined from them depend weakly on the temperature determination. Ions with 
lines in both the infrared and the visible include 0^^, Ne**, S**, and Ar**. 

One important study at optical wavelengths was done by Osterbrock and Flather 
(1959), who compared the intensities of Hg, [0 II] 3726 and 3729 A, and [0 III] 

5007 A out to a distance of 30 arc min from the exciting stars in the Trapezium. 
Spherically symmetric thermal and ionization equilibrium models such as those by 
Rubin (1968, 1985) and Simpson and Rubin (1984) (hereafter SR) do not reproduce 
these observed line ratios; that is, the ratios predicted from the models imply 
too few singly ionized ions in the center and too few doubly ionized ions at 
great distances from the exciting stars. The latter can be explained by the 
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fact that the Orion Nebula is also a reflection nebula, as can be seen from 
continuum observations. In fact, Peimbert (1982) finds that all optical lines 
seen farther than about 12' from the Trapezium are produced entirely by 
scattered light. However, due to a central hole in the dust distribution 
(O'Dell and Hubbard 1965, Simpson 1973, Schiffer and Mathis 197^), the inner 
several arc min of the nebula are negligibly affected by scattering (Peimbert 
1982 ). 

Another reason for studying the Orion Nebula is that the highly obscured, 
compact H II regions that have been modeled by Zeilik (1977), Herter, Heifer, 
and Pipher (1983), Rubin, Hollenbach and Erickson (1983) (hereafter RHE), and SR 
do not show the high ionization that is predicted for the hot massive stars 
required to match the luminosities deduced from radio observations. Possible 
solutions to this discrepancy are that the H II regions are excited by clusters 
of stars (Lacy, Beck, and Geballe 1982, RHE, SR), the ionized gas is mixed with 
dust that softens the stellar spectrum (Herter, Heifer, and Pipher 1983), or 
that the model stellar atmospheres used for the computations do not provide 
realistic representations of the Lyman continuum as a function of frequency and 
stellar effective temperature (RHE and Abbott and Hummer 1985). Because 

the Orion Nebula is nearby and relatively unobscured, we can see the exciting 
stars and directly determine their spectral types. It is excited predominantly 
by the 06p star 6* Ori C (Lee 1968) and several 09.5 stars that provide only a 
small fraction of the total flux. An 06 star supposedly has T^^^-40,000K (Conti 
1973), although more recent work (Pottasch, Wessellus, and van Duinen 1979) 
prefers a lower ^eff ( '38,000K) for this spectral type. Mathis (1985) has 
already suggested that should be lower than 40,000K for Orion, baaed on the 

O^/O*"^ ratio. We know that there is not enough dust to affect the ionizing 
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spectrum in the center where there is the hole in the dust distribution. 
Comparison of new models of the Orion Nebula with the extensive infrared and 
optical observations should elucidate whether the deduced low ionization is a 
general problem, or whether it is confined to obscured dusty H II regions. 

Another problem we wish to address is the abundance of nitrogen relative to 
oxygen. Because nitrogen and oxygen have similar ionization potentials, 

Peimbert and Costero (1969) assumed that their ionization structures are 
identical and that one could correct for the missing N’*"*’ in abundance studies of 
H II regions by multiplying N^/H^ by the ratio (0* + 0**)/ 0*. Peimbert and 
Torres-Peimbert (1977) (hereafter PTP) and Torres-Peimbert, Peimbert, and 
Daltabuit (1980) (hereafter TPPD) studied optical lines throughout the Orion 
Nebula to obtain an N^/0* ratio of 0.10, which they also give as the N/0 ratio. 
However, a line from N** can now be observed in the far infrared at 57pm, as can 
lines of 0** at 52 and 88ym. These three lines are mainly density sensitive 
(the optical lines are more sensitive to temperature), but the density can be 
determined from the [0 III] 52/88 ratio if the lines are optically thin, which 
they are in all but the largest H II regions (Rubin 1968, Simpson 1975). Lester 
et al. (1983) in their study of these three infrared lines in several H II 
regions found that the N**/0** ratio is typically 0.2 or higher, and is 0.30 for 
the radio peak of the Orion Nebula Itself! 

N* does have a lower ionization potential (29eV) than 0* (35eV), and in our 
H II region models (RHE; SR; Rubin 1985) the N"*^* ionization zone is larger than 
the O*"*" zone (except for the dense regions with Tg^^>iio,OOOK). Moreover, in 
regions excited by lower temperature stars and in regions where the matter is 
clumped, the ratio N**/0^* inferred by the technique of Lester et al. (1983) is 
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considerably larger than the input N/0. For 06 stars like 0* Ori C with 
T^^j.-40,000K, the zone is similar to the 0** zone in a nebula of uniform 
density (Rubin 1985). However, the models show that the lower the effective 
temperature of the exciting star, the smaller the 0** region compared to the 
region. Thus high spatial resolution observations should distinguish between an 
enhanced [N III] ratio due to (a) a larger N** zone than the 0 ** zone, 

or (b) some other factor such as different colllsional excitation cross 
sections. 

In this paper, (1) we examine the N**/0** ratio as a function of position 
in the nebula to see if it changes as the models predict. (2) We test the 
previous assumption that the optical line strengths are not affected by 
scattering in the interior several arc minutes by comparing the optical and 
infrared [0 III] and [Ne III] intensities. (3) We investigate the problem of 
the overall degree of ionization being lower than predicted from the spectral 
type of the exciting star observed or inferred from the radio flux - a problem 
that was first discovered in highly obscured compact H II regions. 

An Investigation like this is best carried out by comparing observations 
and predictions of lines from several different elements with significantly 
different ionization potentials (such as nitrogen, oxygen, and neon). To do 
this in a reasonable amount of observing time, we chose not to map the entire 
Orion Nebula, but to observe only one sector. To be sure we had a complete 
sample of the ionization structure as a function of distance from the exciting 
star, we observed in a straight line away from the Trapezium. _We chose the 
strip with the largest number of optical observations so that we could compare 
our detailed theoretical models to the singly ionized elements that have only 
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optical lines. We also wanted to avoid the separate H II regions around the 
e* Ori stars (Peimbert 1982). The final choice was a line due South from the 
Trapezium. This included the ionization front of the bar, which we had to model 
as a shell of higher density, since our models are spherically symmetric. The 
simplification of observing and modeling only one sector has the advantage that 
we are not forced to average over the whole nebula; thus the chemical 
composition and stellar effective temperature herein derived should not be 
affected by the limited solid angle of the investigation. 

The remainder of the paper is divided as follows; In the next section we 
present new observations of the far infrared transitions of 0^*, N**, and Ne** 
in Orion. We then introduce our detailed theoretical models of the nebula in 
Section III. In Section IV the models are compared with our data and other 
published data in a detailed way. Our major conclusions are summarized in 
Section V. 


II. OBSERVATIONS AND DATA ANALYSIS 

The observations were carried out with the 91 cm telescope of the Kuiper 

Airborne Observatory on 1983 November 21 and 1985 February 4, 6, 12, and 14 

using the facility cooled grating spectrometer (COS) described by Erickson et 

al. (1984a, 1985). The aperture profiles, which were measured in the laboratory 

using a small blackbody source, had a FWHM of 37" in 1983 November and 46.5" in 

1985 February. In both cases the beam's effective area corresponded to an 

“8 “8 

equivalent disk -1.04 as large, or 2.7x10 ster in 1983 and 4.3x10 ster in 
1985. The chopper throw was -4' in azimuth (approximately E-W). Sequences of 
four 5 or 10s integrations were taken with the source placed alternately in 
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right and left beam to correct for linear variations in the background flux. In 
1983 November the [0 III] 51.8 and 88.4nm lines and the [N III] 57^3um line were 
observed at the 6 spatial positions (P1-P6) given in Table 1. In 1985 February 
the same three lines were reobserved with improved signal-to-noise at positions 
PI, P4, P5, and P6; the [0 III] 51.8pm line was observed at P7, and the [Ne III] 
36.0pm line was observed at positions P1-P4. The 1985 observations are also 
summarized in Table 1 . 

The positions observed lie in a line starting at 6* Ori C and running 
straight south. Absolute positions were determined to ±5" by offsetting from 
either 6^ Ori C or e* Ori A. The Infrared boresight was verified in flight by 
assuming that the 63.07pm continuum peaks at the Kleinmann-Low nebula (KL) at 
R.A.(1950)-5^32”46.7®, 6(1950)— 5 ®24 '27" (Werner 1982). The total 
bandpass of the six Ge:Ga photoconductor detectors was sufficient to include 
both the line and the adjacent continuum using a single grating position. At 
the different wavelengths observed the spectral resolving power was 3300-4200 in 
1983 November and 2600-3300 in 1985 February. 

On the 1983 November flight the absolute flux calibration and the relative 
detector responses were determined by observing the Mars continuum at each 
wavel^gth. The final reduced spectra are the signals obtained at the map 
positions divided by the signals obtained on the Mars continuum, multiplied by 
the flux from Mars as computed using the model of Simpson et al. (1981), and 
corrected for the difference in the water vapor column density. The 
line-of-slght water vapor was measured to be 16 precipitable microns towards 
Mars and 11 precipitable microns towards KL by observing a water line at 
84.42pm. Since there are no strong features in the vicinity of these lines, the 
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resultant correction to the data was i5%. The [0 III] line fluxes were 
increased by 13 J because of a synchronous demodulation phase adjustment between 
them and the Mara calibration. 

In 1985 February the absolute flux calibration and the relative detector 
responses were determined by observing KL on each flight in the continuum 
adjacent to the lines. Calibration spectra taken before and after the mapping 
observations on a single flight typically agreed to better than ±5?. The change 
in the instrument function between these nearby wavelengths was determined to 
±51t using an onboard calibration source. On 1985 February 4 and 6 the 
line-of-sight water vapor towards KL was measured to be 9 precipitable microns 
by observing water lines at 60.33 and 64.22jim. This water vapor value was 
assumed for 1985 February 12 and 14 as well. The water vapor correction was 
15±8$ for CNe III] 36.0pm, which sits on the shoulder of a relatively strong 
telluric water line, and for the [0 III] and [N III] lines. The final 
reduced spectra are the signals obtained at the map positions divided by the 
corresponding signals obtained on the KL continuum and multiplied by (a) the 
ratio of two water vapor spectra to correct for the difference in transmission, 
(b) the change in the instrument function, and (c) the continuum flux from KL as 
given by Erickson et al. (1981). At both epochs the effect of diffraction in 
the telescope was computed to be ^10% and, hence, was ignored. We estimate the 
overall uncertainty in the absolute calibrations to be ±25>. The agreement 
between the two dates is consistent with this considering the different beam 
sizes. 

The line strengths in Table 1 are the result of a nonlinear least squares 
fit to a flat continuum level and a superimposed Gaussian line profile 
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integrated over a detector bandpass. The four free parameters for the fitted 
curves are the line width, line height, line position, and continuum level. The 
line width provides a measure of the instrumental resolution since the lines are 
unresolved; the line position provides an estimate of the uncertainty in the 
wavelength calibration. For the weaker lines the width and/or position were 
assumed known and the line and continuum strengths were fitted. The error bars 
quoted in Table 1 are one-standard-deviation-of-the-mean statistical errors. 

The additional uncertainty in the ratio of line strengths due to systematic 
effects is estimated to be roughly half the absolute uncertainty quoted above 
for the individual line strengths. 

The [Ne III] 36.0ym line was first detected by Erickson et al. (1984b) in 
M17 and subsequently in the planetary nebula NGC 6543 by Shure et al. (1984). 

The present Orion observations at positions P1-P4 are displayed in Figure 1. 

The line at P2 was observed using two nearby grating positions to obtain a 
better measurement of the continuum. From the present Orion data 3.5km s”^ ) 

and new observations of M 17 (V^gj^»l8km s"^), we find the rest wavelength of the 
line to be 36.013 ± 0.004pm, in agreement with our original detection. 

From the [0 III] 52/88pm line ratio, we have derived the electron densities 
Ng in Table 1 (c. f. Lester et al. 1983) using the collisional excitation cross 
sections of Aggarwal, Baluja, and Tully (1982). Using these densities and the 
57/52um ratio, we derived the values of N**/0** listed in Table 1. With the 
most recent cross sections for 0^* and N** (Nussbaumer and Storey 1979), the 
57/52pm ratio is not very sensitive to density; the error bars.for are 

those of the 57/52uni line strength ratio. As predicted by the models discussed 
below, the ratio does increase slightly with distance from the Trapezium 
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(position PI). The agreement of our data with that of other observers such as 
Lester et al. (1983) and Furniss et al. (1983) is adequate, considering that 
either our spatial positions do not agree or our beam sizes are quite different. 

III. MODELS OF THE ORION NEBULA 

For objects like the Orion Nebula that have a large range of electron 
densities and where different lines of sight sample different ionization zones, 
one can best study the ionization structure and chemical composition by 
comparison with detailed models. Our models do not include far UV or radio 
recombination lines, but focus on the physics of the infrared and optical lines. 

a) Spherical Models 

We first consider traditional spherically symmetric models. Those computed 
here use the H II region model program described by Rubin (1968) and updated by 
Rubin (1983) and SR. The electron densities of the models are those that 
reproduce the radio brightness maps at the same positions as our measurements 
(Table 1 ) and extending further south. The models are intended to reproduce 
only the observations along our cut south from the Trapezium. In Figure 2 we 
have plotted small beam radio brightness temperatures and the Hg measurements of 
PTP scaled to 5 GHz for a line south from the Trapezium. There are two curves 
for Martin and Gull (1976) - the solid line is south from the Trapezium and the 
dashed line is south from the radio peak (<20" west of 0* Ori C). Their beam 
extended 20" in the north-south direction, and shows the bar clearly. A curve 
for the VLA map of Johnston et al. (1983) (M.8 GHz; I6"x13" beam) is also 
plotted. The absolute flux levels are substantially lower, no doubt because the 
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VLA is not sensitive to smooth extended emission. However, the general 
structure of the peak and bar are very similar. The map of Rodriguez and 
Chaisson (1978) (23.** GHz; 80" beam) shows the sharp fall-off beyond the bar but 
none of the structure within the bar. The map of Wilson and Pauls (1984) (also 
23 GHz; 43" beam) is similar. The other maps were made by Johnston and Hobbs 
(1969) (1.6* beam; 9.55mm) and by Schraml and Mezger (1969) (2* beam; 1.95cm). 
These last 4 maps show the outer structure of the H II region that the 
interferometer maps do not, but they are badly degraded by their poor resolution 
within about 4* from the Trapezium. The composite radio map that was used in 
the models reproduces the high radio peak (curve 1 for distances 50.25*), the 
trough and bar of curve 2 out to 2.25*, and then decreases slowly from Tg-60K at 
2.25* to meet the large beam fluxes at 4.5*. (The "trough" is the lower 
intensity region at -I*; the "bar" is at 1. 5-2.0*.) The high density maximum at 
the Trapezium in the model density distribution was used because our 37" and 
46.5" beams include much of the radio peak, even when centered on the Trapezium 
star-0* Ori C. The small disagreement with the Hg measurements of PTP may be 
due to the uncertainties in the extinction correction (approximately one 
magnitude) or to beam size effects. The resultant density distribution adopted 
for the spherically symmetric one component model is shown in Figure 3. which 
includes the effects of cluiq>lng as discussed below. 

The physics of the models is the same as that used by SR with the addition 
of the new collisional excitation cross sections for [Ne II] and [Ne III] 
calculated by Bayes, Saraph, and Seaton (1985) and by Butler and Mendoza (1984), 
respectively, and new charge transfer cross sections for 0^*+He®-*0*+He* by 
Butler, Hell, and Dalgarno (1984). Again we used the line-blanketed LTE model 
atmospheres of Kurucz (1979) for effective temperatures of 40,000K 
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(log g - M.5) and 37"**0,000K (log g ■ 4.0). The luminosities of the stars 
exciting the models were adjusted to give a Strdmgren radius of about 8 arc min, 
which corresponds to the size indicated by the radio maps; the resulting Lyman 

ilO -1 

continuum flux is 10 photons s . The radius of the exciting star for each 

model was adjusted accordingly; they range from 10.2 R for T -«>-40,OOOK to 13.8 R 

© ei r ’ © 

for T^^j.-37,000K, both with log g -4.0. Models with Tg^^-39,OOOK and log g -4.0 
are very similar to the models with T^^j.-40,000K and log g -4.5. Since our 
method determines abundances by comparing model predictions to observations, we 
begin with abundances taken from the literature: He/H-0.10, C/H-3.3x10”**, 
N/H-4.5x10"^, and Ar/H-4.5x1 o“^ (FTP and TPPD) and 0/H-4.0x10~^, Ne/H-8. 1 x1 0"^, 
and S/H»2.2x10~^ (Lester, Dinerstein, and Rank 1979, hereafter LDR). We will 
discuss the uncertainties in these abundances later. For spherical models the 
gas must be clumped, since the observed line ratios that are sensitive to 
density imply much higher densities than implied by the computed ratios for 
models without clumping. We define the clumping factor as the Inverse of the 
filling factor. A clumping factor of 4.0 was derived by considering the ratios 
of the [0 III] 52pm and 88ym lines. It was adjusted until the average of the 
observed ratio equaled the average of the computed ratio for the 6 map 
positions. Figure 4 shows the ionization structure of the spherically symmetric 
model with T^^^-39,OOOK. The effect of the high density bar at 2' on the 
ionization is striking. Moreover, the zone is bigger than the 0^* zone, as 

expected. We note that the fraction of is very small. 

Figure 5 compares observed optical and infrared line intensities and ratios 
plotted against distance from the Trapezium with model predictions. The optical 
line ratios were chosen because they are not as sensitive as the individual line 
intensities to errors in the electron density structure deduced from the 
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deconvolution of the radio maps. We make detailed comparison of 0^ and 
lines because they characterize the ionization structure. We are also 
interested in the ionization structure of nitrogen. Ne** is important because 
its 41 eV ionization potential makes it extremely sensitive to the effective 
temperature of the exciting star. We do not discuss sulfur and argon in this 
paper because we do not have good observational data for all the important 
ionization states. The [S IV] line at 10.5ym has been a notable problem in 
matching theory with observations (RHE and SR) and there are no observations 
with a sufficiently large chopper throw. [Ar II] measurements for the Orion 
Nebula have not been published. Since there are no appropriate hydrogen lines 
to ratio to the infrared lines, the 36pm [Ne III] and the 52ym [0 III] line 
intensities are plotted directly; like the radio continuum, they are direct 
functions of the luminosity of the exciting star as well as its effective 
temperature. 

These spherically symmetric one-component models give a reasonably good 
account of the doubly ionized species (Figures 5a, b, c, d, and e). However, 
they do not describe well the singly ionized species (Figures 5f and g) at 
positions near the center. The observed high intensities imply many more singly 
ionized atoms throughout the H II region than predicted, but particularly in the 
central 2'. Only a model with <37 ,000K could conceivably fit the [0 II] and 

[N II] data. In addition, the density implied by the singly ionized CO II] 
3729/3726 ratio (Figure 5h) of Osterbrock and Flather (1959) and Caplan (1972) 
implies a density roughly 1 .5 times larger than that implied by the model ratio. 
Since most of the 0* is in the low density outer parts, the clumping would have 
to be much larger in the 0^ zone than is required for 0^^ in order to produce 
the observed 3729/3726 ratio. Figures 5i-o will be discussed in the next 
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section, where more detailed comparisons with the observations are made. 


b) Bar/Halo Models 

Some of the shortcomings of the one component models are probably due to 
their geometric simplicity; for example, Zuckerman (1973) suggested that the 
Orion Nebula should not be treated as spherically symmetric. As an alternative 
to spherically symmetric models, we derived a two component bar/halo model of 
the Orion Nebula which we will compare to the standard one component models 
discussed above. The two component model consists of two hemispheres, one in 
front of the star along the line of sight (also called the "halo” model) and one 
in back (the "bar" model) where the erosion of the molecular cloud is occurring. 
Each hemisphere is taken as half of a spherically symmetric model in the 
computation. Figure 6 is a schematic drawing of the bar /halo models. We choose 
this geometry as a first attempt at non-spherlcity because the bar as seen in 
optical photographs is clearly an ionization front (see Goudis 1982 for a 
review) but yet there is low density ionized material extending to great 
distances from the Trapezium. The electron densities are given in Fig. 3. The 
clumping factor derived from the 52/88pm ratio for the bar /halo models is 3.0. 
The front hemisphere has a high density at the center, with density decreasing 
monotonically out to the edge. The back hemisphere starts with the same high 
central density, decreases sharply, but then increases again to an ionization 
front at 0.29pc (2' at SOOpc) from the star. This reproduces the bar in 
projection. The bar/halo geometry reproduces the same radio flux as the one 
component models. The chief advantage of the bar/halo models is that the bar 
hemisphere has its singly ionized ions at high densities concentrated within the 
central 2*. The bar /halo models are also intended to reproduce only the 
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observations along our cut south from the Trapezium. 

Figure 7 shows the predictions of the bar/halo models and the observations. 
(The sharp discontinuities at 2' are due to the edge of the bar.) The bar /halo 
models fit the data for singly ionized species significantly better than the one 
component models. However, it is apparent that even this departure from 
spherical symmetry is not sufficient to explain all aspects of the Orion Nebula; 
more realistic models (which are beyond the scope of this paper) are still 
needed. On the other hand, for doubly ionized species the two types of models 
give very similar answers. Although geometry can certainly explain some of the 
discrepancies of the predicted intensities with the observations, particularly 
the predicted electron density sensitive lines, we show below that other 
discrepancies must have other causes. 

IV. DETAILED COMPARISON OF OBSERVATIONS AND MODELS 

a) Oxygen 

Our data, plotted in Figs. 5a and 7a, show that the bar region is less 
dense than the trough region, because the 52/88um ratio decreases monotonlcally 
with distance from the Trapezium. The bar is less dense in [0 II] as well 
(Simpson 1973). However, Figs. 5b and 7b and Table 1 show that the 52um and 
88ym intensities Increase at the bar, implying a longer pathlength through the 
bar than through the trough, or alternatively, a smaller clumping factor in the 
bar. The clumping factor represents a way of accounting for inhomogeneities in 
models that are unrealistically simple by computational necessity. Our models 
assume that the clumping factor is independent of position. Moreover, because 
of the requirements of spherical symmetry, there is very dense bar material in 
the models at all lines of sight out to 2* (Fig. 6), requiring that the material 
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In the trough line of sight be essentially zero. The lack of agreement between 
the observations and the model predictions in Figs. 5a and 7a shows that 
spherically symmetric geometries or even our bar/halo geometry are not 
appropriate for the Orion Nebula. However, one can partially compensate for the 
inadequate geometry of the models in abundance analyses by averaging nearby 
points. 

The predicted 3729/3726 ratio in the bar /halo models (Figs. 7h) now implies 
too high a density; this is undoubtably the same geometry effect discussed above 
(i.e. the bar is too dense). 

We next compare our infrared data and the optical data of FTP, TPPD, and 

LDR through the ratio [0 III] 52ym/5007A (Figs. 5i and 7i). This ratio was 

constructed from data points that overlapped spatially, albeit imperfectly. The 

ratio is sensitive to electron density N through the 52pm line and to electron 

temperature T^ through the 5007 A line. (Higher density decreases the density 

sensitive line intensities through collisional de-excitation and higher electron 

temperature increases the intensities of~the temperature sensitive forbidden 

lines.) The models show good agreement with the measured line ratios in the 

core, but the model ratios are too high in the halo. This implies that either T^ 

or Ng is too low in the models. The latter is not likely because the bar/halo 

models have higher N^ (Fig. 3) and the agreement is worse. In fact, the 

predicted 52/88pm ratio in the bar /halo models (Fig. 7a) shows that N is too 

high in the halo compared to our observations. The optical ratio 5007/A363A 

(Figs. 5j and 7j) shows that T is too low in all but perhaps the very centers 

© 

of the models. (The 4363A line is more sensitive to electron temperature than 
the 5007A line because it comes from a higher energy level). The bar /halo 
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models are less discrepant in the cores than the one component models. Most of 
the discrepancy is probably due to [0 III] 4363A (a favorite electron 
temperature indicator), because the 52/5007 and 5007/HB ratios, and the log of 
the 52pm intensities (Figs. 5i and 7i, 5d and 7d, and 5b and 7b, respectively) 
all agree better with the observations than does the 5007/**363 ratio. If the 
problem is the electron temperature in the models (rather than the cross section 
for the 1I363A line), would have to be raised from ~7000K at 2.5' from the 
center to over 8000K to provide the required agreement (but see the later 
section on neon). 

The optical ratio [0 II] (3726+3729 )/Hg (Figs. 5f and 7f) shows that there 
is not enough singly ionized oxygen, particularly in the one component model. 

The bar /halo model gives a better fit but not in the central 1 ' . The inadequacy 
of the model predictions cannot entirely be attributed to abundance, because 
plots of the [0 III] 5007/Hg (Figs. 5d and 7d) and the log of the 52um intensity 
iFlgs. 5b and 7b) show better agreement with the observations. Some of the 
excess [0 II] 3727 A observed in the center may be due to scattering; however, 
this must be a small effect because the absolute flux at 3727 A is so much larger 
in the center than at 4' from the center. A more likely explanation is that the 
assumed stellar atmospheres are inappropriate. An exciting star with Tg^j.<37,000K 
seems to be necessary for agreement with 0^, although 0^^ agreement is better 
with T^j.j.» 37-38 ,000K. If the exciting star in the models were less luminous 
(i.e. had a smaller radius), the 0** Strdmgren spheres would be smaller (too 
small for the radio maps) and the strong decrease in [0 III] Intensities seen in 
Figs. 5b, 5d, 7b, and 7d that marks the StrOmgren sphere would come at smaller 
distances from the Trapezium. This would help the 0*^ agreement, but not the 0^ 
agreement. 
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b) Nitrogen 


As with singly ionized oxygen, the models need more singly ionized 
nitrogen, as can be seen in the 6584/H6 ratio in Figs. 5g and 7g when the 
stellar effective temperature is >37»OOOK. Although scattering should not 
contribute as much to the intensity of the 6584A line as to that of the 3727A 
line (because the cross section increases as 1/A), the 6584/Hg excess appears to 
be as large as the (3726+3729 )/H0 excess (Figs. 5f and 7f). The electron 
temperature, as determined from the 6584/5755 ratio (Figs. 5k and 7k), could be 
higher for the one component model, but is satisfactory for the bar /halo models, 
particularly in the core. Even if the electron temperatures of the models were 
increased within reasonable limits, thereby increasing the [N II] 6584/Hg ratio 
(Figs. 5g and 7g), the models would still require more N*. Since the input N/0 
ratio (0.1125) is derived from N'*^/0*, the agreement in Figs. 5i 2 ind 7H is 
expected to be good; in fact the ratio could be 15$ larger in the bar /halo 
model, and 35$ higher in the one component model. The N/0 ratio for the sun is 
0.1175 (Aller 1984). 

On the other hand, the observed 57/52um ratio (Figs. 5m and 7m) of N*'*'/0'^'*’ 
is about a factor of 2 larger than any model ratio in spite of the fact that the 
N^^ zone is bigger than the 0 ** zone for all the models. Likewise, the derived 
N^’*’/0'*’^ ratio is larger than the optical N*/0* ratio, a problem common to other 
H II regions (Lester et al. 1983, Dinerstein et al. 1984). Previous proposed 
solutions for this discrepancy Included the effects on the ionization 
equilibrium of having a low effective temperature for the exciting star (Rubin 
1983, Dinerstein et al. 1984), a large amount of collisional de-excitation in 
the dense gas near the star where the ionization levels are similar (this gives 
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more weight to the low density material far from the star where nitrogen is more 
ionized than oxygen), or optical depth effects in the infrared lines (Lester et 
al. 1983). The models show that the differences in the N** and 0*^ sampled at 
various lines of sight appear only at large distances (impact parameters) from 
the center (>2' to 3*). This eliminates the first two proposed solutions. 

Moreover, the agreement of the 52pm fluxes with the models (Figs. 5b and 7b) and 
with the 5007 A line (Figs. 5i and 7i) shows that the optical depth at 52pm (and 
88pm) is small. We clearly see the advantage of mapping the nebula, because we 
can immediately eliminate all explanations due to ionization equilibrium! 

Surviving potential solutions to this problem include errors in the 
electron density or in the atomic data. Since [0 III] 52pm and [N III] 57pm 
have almost the same sensitivity to electron density when the current cross 
sections are used, errors in the assumed electron densities are probably not the 
solution. If the cross sections are in error, more than one cross section must 
be wrong. If just the [N III] 57pm cross section were too low, raising it by a 
factor of two would change the intensity by only a small amount because the 
densities in the models are already so high that the 57pm line is collisionally 
de-excited to a large extent. If the cross section for the [0 III] 52pm line 
were decreased by a factor of two, the 5007A cross section would also have to be 
correspondingly decreased because the 52pm/5007A ratio agrees with the 
observations (Fig. 7i). The agreement with the 5007/4363 ratio (Fig. 7J) would 
then be improved and the derived oxygen abundance would be closer to that of the 
sun (0/H-8.1X10"'*, Aller 1984). 

It seems most probable that the nitrogen abundance really should be a factor of 
two larger than that used in our models, in spite of the resulting disagreement 
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with the optically derived N*/0* ratio (and the solar N/0 ratio). If the [N II] 
cross sections need to be revised, those for [0 II] 3726A and 3729A would have 
to be increased as well because of the apparent agreement shown in Figs. 5 % and 
74. This would give better agreement with the [0 II] 3727/HB ratios (Figs. 5f 
and 7f) as well as the [N II] 6584/HB ratios (Figs. 5g and 7g). 

c) Neon 

The [Ne III] 3869/Hg ratio (Figs. 5e and 7e) and the plot of the log of the 
36um intensity (Figs. 5c and 7c) both require Tgj.j.«39-40,000K. One could say 
that the [0 III] 52pm line could fit a higher temperature star too if only the 
inner 4 data points (to 3') were considered. However, there is little 
difference between the 5 models for the 52pm line in this region, and only the 
lower T^^j. is consistent with the outer data points. On the other hand, there 
is much greater discrimination between the various models in the inner 3' region 
for the [Ne III] 36pm intensities than for those of [0 III] 52pm. Clearly, the 
[Ne III] 36pm line cannot be fitted well with the 36-38, OOOK models that are 
required for oxygen and nitrogen. The optical [Ne III] and [0 III] observations 
give the same results as the infrared, thus showing that the problem is the 
ionization equilibrium (or stellar atmosphere), and probably not the collisional 
excitation cross sections. However, the latter may be responsible for the 
disagreement in the electron temperatures derived from each element (the oxygen 
and nitrogen lines showed that the electron temperatures are too low in the 
models, but the [Ne III] 36pm/3869A ratio (Figs. 5n and 7n) shows that the 
electron temperature is too high in the models). 

Unlike nitrogen and oxygen which are mostly doubly ionized, most of the 
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neon in the nndels is in the Ne* ionization state. The only forbidden line of Ne* 
is at 12.8)iin, where most of the ground based telescopes chop with a throw 
smaller than the nebula. The [Ne II] observation with the largest chopper throw 
(100”) is that of LfCSter ( 1 982;pers. comm). He measured fluxes of 8~12xl0 cm ^ in 5 
positions on a grid just west of 0* Ori C with a 10" beam. His [Ne II] fluxes 
are plotted in Figs. 5o and 7o. The Ne^ intensities are well fitted by models 
with Tg^j.-37-38,OOOK. As with and 0*, there is too much singly ionized neon 
for the higher models implied by the doubly ionized species. For neon, a 

compromise with T^^^-39,OOOK and the Ne/H abundance higher by 20-50% fits both 
ionization states reasonably well. However, this does not agree with the 

oxygen and nitrogen data. A better fit might be obtained if the emergent 
stellar flux resembled that of the Kurucz 37,OOOK, log g - ^.0 model (Kurucz 
1981) from 13.6eV to 41 eV, but with the flux decreasing less rapidly from 41 eV 
to 54eV. 


V. SUMMARY 

We present new observations of the Infrared lines of [N III] 57ym and 
[0 III] 52ym and 88ym at 6 positions and [Ne III] 36ym at 4 positions in a line 
south from the Trapezium in the Orion Nebula. From the data we derive electron 
densities and N^^/0^^ ratios at all 6 positions. 

We also present new models of the Orion Nebula, both spherically symmetric 
and two component bar /halo models, in which the electron density distributions 
are constrained to reproduce the high resolution radio surface brightnesses at 
the same positions as our data. The stellar effective temperatures range from 
37.000K to 40,000K with log g-4.0 or 4.5. The stellar Lyman continuum fluxes 

i|Q .1 

are 10 photons s . The two-component bar /halo models give a somewhat better 
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account of the singly Ionized atoms. 


We compare the predicted intensities from the models to infrared and 
optical data to investigate a number of issues: 

a) Elemental Abundances 

The abundance ratio 0/H»4.0xl0 used in the models, which is approximately 
half solar, agrees satisfactorily with the observations. However, given the 
current cross sections, there is enough discrepancy between the 0* and the 0^* 
ionization structure that the uncertainty is large. The Ne/H abundance ratio is 
almost certainly larger than the 8.1x10”^ used in the models; 1.3x10~** is 
probably a better value, although then the Ne/0 ratio would be larger than the 
ratio found in the sun and most planetary nebulae (Aller and Czyzak 1983, Aller 
1984). Based on our [N III] and [0 III] measurements we derive a N/0 ratio of 
0.21. The ratio from the optical [N II] and [0 II] lines is N/0-0.11. [N III] 
and [0 III] should be a better indicator of N/0 because 0* Ori C is hot enough 
that nitrogen and oxygen are mostly doubly ionized in the nebula. Because of 
this, we revise the nitrogen abundance upward to N/H-8.4x10~^. Thus we find 
that nitrogen and neon are approximately solar in abundance, but compared to the 
sun, oxygen is underabundant. It is conceivable that the missing oxygen is tied 
up in dust grains, although this hypothesis has previously been rejected for the 
Orion Nebula (e.g. Peimbert 1982, Mathis 1985). 

b) Stellar Effective Temperature 

The [0 III] and [N III] lines imply that the effective temperature of the 
exciting star is cooler than 37,000K. However, the [Ne III] lines, which are 
more sensitive to T^^^, lead us to conclude that the effective temperature is 


-23 


S39,OOOK for either model geometry. With no changes in abundances, the observed 
[Ne II], [N II], and [0 II] line fluxes would require that be 37,000K, 
<37,OOOK, and <37,OOOK, respectively. However, with the higher abundances for 
neon and nitrogen discussed in the previous section, the singly and doubly 
ionized species can both be fitted by temperatures of 37,OOOK for nitrogen and 
39,OOOK for neon. Oxygen cannot be fitted easily, but certainly implies 
Teff i 37,OOOK. Thus in this relatively unobscured H II region, the 
derived from nitrogen and oxygen is significantly lower than the derived 

from the neon ionization. As mentioned in the introduction, in obscured H II 
regions the stellar effective temperature derived from the sulfur and argon 
ionization equilibrium is much lower than that estimated from the radio 
luminosity. From our additional data on neon and oir mapping of the Orion 
Nebula in oxygen and nitrogen, we conclude that the presence of dust does not 
account for the inconsistencies in the methods of derivation of stellar 
effective temperatures. In the future, we plan to investigate other stellar 
spectra in the manner of RHE, such as that of a 37,OOOK star but with additional 
flux beyond 41 eV to produce additional Ne*"*^. In this way, we hope to find a 
model that better fits the observed spectrum of the Orion Nebula. 

c) Scattering by Dust 

All the models show that the ratio of the singly ionized forbidden lines to 
Hg increases faster than observed with distance from the center. If some of the 
visual flux observed in the center of the Orion Nebula is due to light scattered 
into the beam from positions far 02') from the center, there would be less need 
for a low ^eff to produce the large number of singly ionized atoms observed. 
Scattering should be much more important for [0 II] 3727A than for [N II] 6584A. 
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However, the two ions show similar discrepancies, leading us to conclude that 
scattering is not the chief cause of the large singly ionized line flux within 
2’ of the center, but rather, it is the ionization equilibria and abundances. 
However, scattering could still affect the slope of the line/HB ratio at the 
center. On the other hand, our infrared oxygen and neon data show good 
agreement with optical data and demonstrate that scattering is also not 
important for doubly ionized ions in the central 4' of the nebula. 

d) Electron Densities 

The measured [0 III] 52/88 ym line ratios show that the enhanced 
intensities at the ionization bar are due to geometry effects rather than 
increased electron densities. We therefore do not expect our present models to 
reproduce the details of the density-sensitive line ratios in the vicinity of 
the bar. An average clumping factor for the models was derived from the far 
infrared [0 III] lines, and the density variations in the models ^ reproduce 
the general Increase in the line ratios of singly ionized ions (e.g. [0 II] 
3729 / 3726 ) with distance. The detailed agreement might be improved by using a 
different geometry /density distribution such as a blister model (Rubin 1984). 

We note however that the blister geometry will not account for the enhanced 
intensities near the bar. 


e) Electron Temperature 

The model electron temperatures are too high in the very center (Ne*^ 
zone), but low in the 0^^ and N* zones. The Ne** and 0^"^. zones, of course, 
overlap. For 0^*, the model temperatures are worse for [0 III] 5007/4363 than 
for [0 III] 52ym/5007A. Electron temperatures are not very sensitive to T^^^ 


-25 



for the range of considered here. These discrepencies are possibly due to 
errors in the atomic parameters. 
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INFRARED LINE OBSERVATIONS OF THE ORION NEBULA 
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FIGURE CAPTIONS 


Fig. 1. The [Ne III] lines at 36ym at each of the four positions P1-P4 (top to 
bottom). 

Fig. 2. The radio brightness temperature scaled to 5 GHz for a cut south from 
the Trapezium. The H8 measurements of PTP as scaled to 5 GHz are also 
plotted. 

Fig. 3. The electron densities in the clumps as a function of distance from the 
exciting star. The units of the abcissa are 0.l45pc, which corresponds to 
1 ' at a distance of 500pc. 

Fig. 4. The ionization structure in nitrogen, oxygen and neon as a function of 
distance from the center for the one component model with T^^^-39,OOOK. 

The units of the abcissa are 0.l45pc, which corresponds to 1 • at a distance 
of 500pc. 

Fig. 5. Line intensities, or intensity ratios, for the five different models as 
functions of distance from the center of the model. The observed data are 
also plotted as functions of distance from the Trapezium. With exceptions 
as noted, the squares are our data of November, 1983» and the triangles are 
our data of February, 1985. Where no error bars are plotted for our data, 
the error bars are smaller than the size of the plotted squares or 
triangles. The circles are the visible data of PTP and TPPD, and the 
diamonds are the visible data of LDR. Fig. 5a. The ratio of [0 III] 
52/88um. 

Fig. 5b. The log of the [0 III] 52ym Intensities. 

Fig. 5c. The log of the [Ne III] 36ym intensities. 

Fig. 5d. The ratio of [0 III] 5007A/H8. 
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Fig. 5e. The ratio of [Ne III] 3869A/He. 

Fig. 5f. The ratio of [0 II] (3726+3729 )A/Hb. 

Fig. 5g. The ratio of [N II] 6584A/H0. The triangle is from O'Dell and 

Hubbard (1965). 

Fig. 5h. The ratio of [0 II] 3729/3726A. The circles are from Osterbrock and 
Flather (1959), the squares are from Simpson (1973), and the x's are from 
Caplan (1972). 

Fig. 5i. The ratio of [0 III] 52um/5007A. 

Fig. 5j. The ratio of [0 III] 5007/^363A. 

Fig. 5k. The ratio of [N II] 6584/5755A. 

Fig. 51 . The ratio of [N II] 6584A/CO II] (3726+3729 )A. 

Fig. 5m. The ratio of [N III] 57um/[0 III] 52pm. 

Fig. 5n. The ratio of [Ne III] 36pm/3869A. 

Fig. 5o. The log of the Intensity of [Ne II] 12.8ym. 

Fig. 6. Schematic drawing of the bar /halo models. 

Fig. 7. The same as Fig. 5 for the bar/halo models. Fig. 7a. The ratio of 
[0 III] 52/88pm. 

Fig. 7b. The log of the [0 III] 52pm intensities. 

Fig. 7c. The log of the [Ne III] 36pm intensities. 

Fig. 7d. The ratio of [0 III] 5007A/H6. 

Fig. 7e. The ratio of [Ne III] 3869A/H0. 

Fig. 7f. The ratio of [0 II] (3726+3729)A/HB. 

Fig. 7g. The ratio of [N II] 6584A/H0. 

Fig. 7h. The ratio of [0 II] 3729/3726A. 

Fig. 7i. The ratio of [0 III] 52pm/5007A. 

Fig. 7j. The ratio of [0 III] 5007/4363A. 
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Fig. 7k. The ratio of [N II] 6584/5755A. 

Fig. 1 %. The ratio of [N II] 658i»A/[0 II] (3726+3729)A. 

Fig. 7m. The ratio of [N III] 57pm/[0 III] 52ym. 

Fig. 7n. The ratio of [Ne III] 36ym/3869A. 

Fig. 7o. The log of the intensity of [Ne II] 12.8ym. 
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